The thermodynamic spin magnetization of strongly correlated 2D electrons in a 

silicon inversion layer 
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A novel method, invented to measure the minute thermodynamic magnetization of dilute two 
dimensional fermions, is applied to electrons in a silicon inversion layer. The interplay between the 
ferromagnetic interaction and disorder enhances the low temperature susceptibility up to 7.5 folds 
compared with the Pauli susceptibility of non-interacting electrons. The magnetization peaks in 
the vicinity of the density, where transition to strong localization takes place. At the same density, 
the susceptibility approaches the free spins value (Curie susceptibility) , indicating an almost perfect 
compensation of the kinetic energy toll associated with spin polarization by the energy gained 
from the Coulomb correlation. Yet, the balance favors a paramagnetic phase over spontaneous 
magnetization in the whole density range. 



INTRODUCTION 



The nature of the ground state of degenerate two di- 
mensional (2D) fermions at zero magnetic field is an out- 
standing open problem, which has not been deciphered 
despite decades of research. In the absence of disorder the 
ground state is believed to be determined by an interplay 
between the kinetic energy, Ep, and the inter-particle 
interaction energy, Ec = e^/na, where a = (tto)"^/^ is 
the inter particle distance, n is the areal particle den- 
sity, and K is the host dielectric constant. The relative 
importance of the two energy scales is characterized by 
Ts = a/oo, with flo being the Bohr radius. For electrons 
in a single band ~ Ec/Ep, while for the (100) sur- 
face of silicon rg — Ec/lEp due to the two- fold valley 
degeneracy. At very high densities (rj <C 1) a 2D sys- 
tem approaches the non-interacting degenerate gas para- 
magnetic limit, characterized by the Pauli susceptibility 
Xo- As the density is reduced, the growing ferromagnetic 
correlations lead to substantial enhancement of the spin 
susceptibility x- The system is predicted to remain para- 
magnetic up to « 20 25, where numerical calcula- 
tions find a quantum phase transition to a ferromag- 
netic liquid phase 0|. At lower density, fa 37 the 
Coulomb correlations are predicted to lead via another 
phase transition to a quantum Wigner crystal with frus- 
trated antiferromagnetic spin arrangement followed 
by transition to a ferromagnetic arrangement at an even 
lower density The energy balance between the ferro 
and paramagnetic states is very subtle and the density 
window where ferromagnetism may take place is small 
P, [^. Such a ferromagnetic phase has never been ob- 
served experimentally. The situation is fundamentally 
complicated by the unavoidable disorder present in any 
realistic system. In the absence of Coulomb interactions 
all wave functions of a 2D system are believed to be expo- 
nentially localized J_] . Localization modifies the Coulomb 
interaction dramatically in the low density limit. The 
interplay between kinetic energy, interaction, and disor- 



der was worked out theoretically for the case of relatively 
weak disorder 0,0,^3- '^^^ found that the interaction 
suppresses the localizing effect of disorder, especially in 
the presence of valley degeneracy fll| . Yet, at low enough 
densities disorder prevails and localization always com- 
mences. 

Notwithstanding the substantial research done thus 
far, there is presently no agreed picture of the phase dia- 
gram corresponding to a realistic 2D fermion system. It 
is clear that the spin degree of freedom plays a crucial 
role in the low density regime, n < 2 x 10^^ cm^^, but 
the minute total magnetic moment pertaining to such 
a small number of spins has hindered, thus far, any di- 
rect measurement of the thermodynamic spin magneti- 
zation. Present estimates of the 2DEG magnetization in 
silicon rely on susceptibility data obtained from trans- 
port measurements, either Shubnikov-de Haas (Sh-dH) 
oscillations in a tilted magnetic field ll^l 0| 

or saturation of the magnetoresistance in an in-plane 
field 0, 0| . The two approaches led to contradicting 
conclusions. While the magnetoresistance data were in- 
terpreted as indicating the long awaited Bloch-Stoner Q 
instability at the critical density for the metal-insulator 
transition, analysis of Sh-dH oscillations points against 
such instability 18]. 

At the heart of the present manuscript is a novel 
method invented to measure the thermodynamic magne- 
tization directly. We apply the method to a high mobility 
2D electron layer in silicon. In particular, we find that, as 
the density is reduced, the weak field spin susceptibility 
is progressively enhanced up to 7.5 xoj but the ferromag- 
netic instability is never realized. The system turns in- 
sulating before it polarizes and electron localization then 
leads to a reduction in the Coulomb interaction. The 
localization transition is thus characterized by a sharp 
cusp in magnetization. Interestingly, we find indications 
for localized magnetic moments in coexistence with the 
itinerant electrons, even at high carrier densities. 
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METHOD, SAMPLES, AND EXPERIMENTAL 
SETUP 



The experimental setup is presented in Fig. ^ An ex- 
ternal bias, Vg, sets a constant electrochemical potential 
difference between the gate and the 2D channel equal to 
the sum of the electrostatic potential difference, (j), and 
the difference between the aluminum gate and the 2DEG 
work- functions, Wai and W2D, respectively 



eVa = e^{n) + Wai - W2D{n, B). 



(1) 



Modulation of the in-plane magnetic field by an auxiliary 
coil at a frequency w modulates the chemical potential of 
the 2DEG and, hence, W2D {Wai modulation is negligi- 
ble). Since Vq is kept constant, the differential of Eq. 
vanishes. The 2DEG chemical potential, /i, equals the 
Si — Si02 band discontinuity minus W2D (Fig-Dl- Con- 
sequently, one obtains 
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where {ed4)/ dn + dji/ dn) / = C^^ is the independently 
measured inverse capacitance per unit area, compris- 
ing the geometrical and chemical potential contributions. 
The latter contribution includes well width and interac- 
tion effects. 

In terms of the induced current, 51 , and the magnetic 
field modulation, 5B, one obtains 



dfi ieSI 
dB ^ ^ CudB " 



(4) 



Since the 2D layer thickness and the screening length 
are minuscule compared with the oxide thickness, the 
capacitance is close to the geometrical one, and hence, 
constant to within 1% in the whole density range. Us- 
ing one of Maxwell's relations, dM/dn = —d^/dB, we 
obtain dM/dn and integrate it numerically with respect 
to n to derive the magnetization M{B, n). The magnetic 
susceptibility x is calculated from the slope of M{B,n) 
versus B at small fields. An additional constant field, 
induced by the main coil, facilitates magnetization mea- 
surements at finite magnetic field. 

While the method is conceptually straightforward, its 
realization is demanding since the current induced by 
the field modulation is typically on the order of 10~^^ A, 
while the spurious current induced in all wire loops by 
the ac magnetic field, and even more so by mechanical 
vibrations of the sample in the dc magnetic field, are 
potentially larger by several orders of magnitude. The 
induced currents were minimized in the experiment by 
careful compensation of all loops. Mechanical vibrations 
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FIG. 1: Magnetization measurement setup and band diagram 
of the 2D confining potential 



were minimized by rebuilding the relevant parts of the 
refrigerator to achieve high enough mechanical rigidity. 
After building several prototypes we were able to drive all 
mechanical resonances to frequencies considerably higher 
than the field modulation frequency, and hence, eliminate 
the mechanical vibrations at the measurement frequency. 

The same setup can be used to measure the much 
larger orbital magnetization in a perpendicular magnetic 
field. The method's sensitivity scales with sample's area 
and magnetic field modulation amplitude. For the 4 mm^ 
sample and 0.03 T rms field modulation used here, it 
was about 10~^^ J/T, comparable to the best sensitiv- 
ity achieved with SQUID-based magnetometers (a 
review of the current state of the art in magnetization 
measurements can be found in |20() . The advantage of 
our (and'l?) method is in its applicability to arbitrary 
magnetic fields and temperatures, as well as to a wide 
range of conductivities. The extraction of the magnetiza- 
tion from Sh-dH oscillations, on the other hand, requires 
perpendicular magnetic fields, low temperatures and high 
enough mobilities. The most interesting regime, the tran- 
sition to strong localization is, hence, at the limit of its 
reach. 

The samples used in the experiment were similar to 
those used in Refs. 0, They consisted of 5 mm long 
0.8 mm wide Hall bars with 2.5 mm separation between 
the potential probes. The oxide was 200 nm thick, lead- 
ing to C = 678 pF device capacitance. We applied —15 V 
bias to the substrate in order to minimize the contacts 
resistance 23]. The maximal mobility under this bias 
reached 17, OOOcm^/Vs . An alternating magnetic field 
of typically 100-300 G (rms) and / = w/27r = 5 20 Hz 
was applied parallel to the layer along with a desired 
constant field. A preamplifier with ^ 2 fA/^/Hz cur- 
rent noise was used to measure the current and bias the 
gate (Fig. 1). 

Unlike Sh-dH based measurements, our method is sen- 
sitive to the total thermodynamic magnetization com- 
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prising the spin part as well as the diamagnetic orbital 
contribution due to the finite ('^ 50 A) thickness of the 
2D layer (Fig. . Localized states also contribute to the 
measured magnetization, as long as they exchange parti- 
cles with the 2DEG at a rate faster than oj. 



EXPERIMENTAL RESULTS AND DISCUSSION 

The measured dM/dn at 9 T magnetic field and T = 
100 mK is depicted by dots in Fig. Et- The smooth 
solid line depicts the same quantity as extracted from 
Sh-dH data [l5|. The difference between the two curves 
is attributed to the diamagnetic shift due to the sub- 
band energy level (eq in Fig. ^| dependence upon the 
magnetic field and to the presence of localized spins. 
Both effects do not appear in Sh-dH oscillations. While 
the localized spins are certainly relevant to the study 
of spin magnetization in real samples, the diamagnetic 
part reflects an orbital effect, which is outside the scope 
of our interest. To estimate the latter contribution we 
assume that the Sh-dH measurements at high densities 
(say n > 5 X 10"'^^ cm~^), where the number of local- 
ized spins is small, give the spin magnetization correctly. 
The diamagnetic contribution is then given by the differ- 
ence between our measured thermodynamic magnetiza- 
tion and the one extracted from the Sh-dH data. At zero 
density, on the other hand, one can calculate the dia- 
magnetic shift in a single-particle picture. To complete 
the estimate at intermediate densities we interpolate be- 
tween the two limits to obtain the dashed line in Fig. 
The spin magnetization in the whole density range is 
obtained by subtraction of the diamagnetic contribution 
(dashed line) from the measured data (dots). The effect 
of the magnetic field, even at 9 T, on the subband en- 
ergy Eq is much smaller than the inter-subband spacing. 
Therefore, the diamagnetic contribution to the magne- 
tization should depend linearly upon magnetic field, in 
accordance with our high density data. The overall dia- 
magnetic contribution in the low-density range is small 
compared with the spin contribution. Moreover, it varies 
slowly with density. The extracted spin magnetization is 
therefore only slightly affected by the details of the inter- 
polation procedure. Yet, for n < 2 x 10^^ cm~^ we find 
that the saturation value of the extracted spin magneti- 
zation dM/dn (solid line in Fig.|2K) is lower by « 10 % 
than the one Bohr magneton per electron, expected for 
full polarization. Since all spins are likely to be polar- 
ized at low density and 9 T, we attribute the discrepancy 
to an underestimate of the diamagnetic contribution at 
low densities. The error in our measured data is much 
smaller than 10 %. The magnetization values presented 
below may, hence, underestimate the actual magnetiza- 
tion at low densities by up to 2 x 10^ /is cm~^ per Tesla. 
This uncertainty is immaterial for our conclusions. The 
rest of the paper focuses on the spin contribution ob- 
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FIG. 2: (a) Total and spin dM/dn. Dots - measured total 
dM/dn, dashed line - diamagnetic contribution ^ th ick smooth 
line - spin dM/dn extracted from Sh-dH data [l^. thin solid 
line - spin dM/dn extracted by subtraction of the diamag- 
netic contribution from the total dM/dn. (b) Spin magneti- 
zation obtained by integration of dM/dn data. The dashed 
line demonstrates that the curves extrapolate to zero M at 
vanishing density, as they should. 



tained by the above procedure. 

The spin magnetization (Fig. \Bp) at a given field is ob- 
tained by numerical integration of the extracted dM/dn 
values with respect to n. Since the magnetization can 
be measured only above a certain density, for which the 
sample resistance is lower than ~ 1 MSI, the integration 
cannot start from zero density, where M — 0. Conse- 
quently, our integration yields the magnetization up to 
a constant, which is chosen so that the magnetization at 
high densities equals the values extracted from the Sh- 
dH oscillations. We neglect the small number of localized 
spins, which are present even at high densities. The fact 
that the resulting curves at all temperatures extrapolate 
to practically zero magnetization at n = (dashed lines 
in Fig. 13)) confirms that the integration constants are 
chosen properly. 

The spin magnetization at various magnetic fields and 
four temperatures is depicted in Fig. For each mag- 
netic field the curves with higher magnetization values 
correspond to lower temperatures. The thick blue line 
corresponds to full polarization of all carriers at a given 
density and the thick red line to the zero temperature 
magnetization of a non-interacting degenerate electron 
gas at i? = 6 T. The empty circles denote for each mag- 
netic field the critical density nc{B), which separates the 
metallic regime from the insulating one (2^. At higher 
densities the resistance decreases as the temperature is 
reduced (metallic behavior) while at lower densities it in- 
creases (insulator). Whether the metallic behavior indi- 
cates a true 2D metal or merely finite temperature trans- 
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FIG. 4: Inverse susceptibility as determined from M{B) at 
B — 0.7 T. Experimental points from bottom to top corre- 
spond to densities 0.8-^6 x lO^^cm"^ in 4 x 10^° cm~^ steps. 
The thick straight line depicts Curie law and the dashed line 
marks T = (g/iB/kB) x 0.7 T. The experimental points at 
n = 8 X 10^" are connected by a thick line for comparison 
with the expectation for non-interacting electrons of the same 
density. 



FIG. 3: (a) Spin magnetization as a function of density at 
different magnetic fields and temperatures 0.2, 0.8, 2.5 and 
4.2 K; higher magnetization corresponds to lower tempera- 
ture. Critical densities, ric, are marked by circles. Thick blue 
line - full magnetization, thick red line - magnetization of a 
degenerate ideal electron gas at B = 6 T. (b) Maximal spin 
magnetization and spin magnetization at the critical densities 
plotted against magnetic field. Dashed line - extrapolation 
from high magnetic fields. 



port through localized states with long enough localiza- 
tion length is presently an open question. It is clear, 
though, that the insulating regime corresponds to local- 
ized states (either in the sense of percolation or in the 
sense of exponentially decaying wave functions) with pro- 
gressively smaller localization lengths at lower densities. 
At high magnetic fields, full spin alignment persists up 
to densities considerably higher than those predicted for 
non-interacting electrons (compare the non-interacting 
and the experimental curves for _B = 6 T). Curiously, for 
all magnetic fields the magnetization reaches its maxi- 
mal value at densities only slightly lower than the criti- 
cal density ric- As more carriers are added to the layer, 
the total magnetization is monotonically reduced. The 
large negative slope of the curves in the vicinity of Uc 
indicates that the added delocalized electrons prefer to 
occupy the upper spin subband. At still higher densities 
the magnetization is further reduced towards the respec- 
tive non-interacting values. 

Fig. depicts the maximal magnetization as well as 
the magnetization at the critical densities versus mag- 



netic field. The data set an upper limit on the zero field 
polarization at the critical density, ric ~ 1.25x 10^^ cm~^, 
to less than 2 x 10^° /ib cm~^. Our data, hence, point 
against Stoner instability in our samples. 

Within effective medium theory the inverse suscep- 
tibility per particle assumes Curie- Weiss form, x = 
fi^/kB{T — Tc). The value of Tc in this approximation 
provides an intuitive measure for the combined effects of 
the kinetic energy and interaction. In particular, a ferro- 
magnetic instability requires positive Tc. The paramag- 
netic nature of the 2D layer should, hence, be reflected in 
the temperature dependence of the susceptibility x- The 
inverse susceptibility, determined from U = 0.7 T spin 
magnetization, normalized by the electron density and 
expressed in Kelvins is depicted as a function of temper- 
ature in Fig. 2] For all densities the inverse susceptibil- 
ity per spin is larger (negative Tc) than the Curie value, 
X~^ = kBT/^^, indicating that in the balance between 
the Coulomb energy gain and the kinetic energy toll as- 
sociated with spin polarization in the system, the latter 
wins. Yet, at the lowest densities the victory is marginal, 
Tc = 0.2K. To appreciate the almost perfect balancing 
of the kinetic energy by the interaction, we compare the 
data for n = 8 X 10^° cm~^ (thick line connecting the 
data points in Fig. 0J) to the theoretical inverse suscep- 
tibility of a non-interacting 2D Fermi gas of the same 
density. The difference between the two curves reflects 
the effects of the ferromagnetic interaction and disorder, 
which are absent in an ideal non-interacting gas. Re- 
markably, the susceptibility measured at densities just 
below ric, approaches the free spin one (Curie law) very 
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FIG. 5: Magnetization as a function of magnetic field at T = 
100 mK and T = 2.5 K. Densities are given in lO^'^ cm~^ units. 
Bold line - Pauli magnetization for non-interacting degenerate 
fermions. As the temperature and magnetic field are reduced, 
the magnetization becomes increasingly nonlinear. 



sufficiently low magnetic fields are required (less than 
0.7 T in our case) to determine the "zero" field suscep- 
tibihty. At T = 100 mK and n = 1.25 x 10" cni'^ the 
strong Coulomb interaction is manifested in a 7.5 folds 
enhancement of the susceptibility compared with non- 
interacting electrons. This susceptibility is twice as large 
as the value extracted from the Sh-dH data for the 
same density. We attribute the difference to the localized 
states which persist into the metallic phase, n > ric, and 
are not sensed by the Sh-dH analysis. The weak inter- 
action between such spins should result in a very large 
susceptibility at low temperatures. Indeed, as evident 
from Figs. 0] and El the weak field susceptibility is sensi- 
tive to temperature down to 100 mK. For stronger fields, 
for which the Zeeman energy exceeds the temperature, 
almost all spins are polarized and susceptibility depends 
very weakly on temperature. 



COMPARISON WITH RECENT SEARCHES FOR 
THE STONER INSTABILITY 



closely, implying that the kinetic energy is almost per- 
fectly compensated by the interaction. Yet, the former 
wins and paramagnetism prevails. Since we believe that 
the free spin-like susceptibility near the critical density is 
generic, rather than fortuitous, we propose that the lo- 
calization transition is driven, in addition to disorder, by 
the strong exchange interaction, which promotes localiza- 
tion through the Pauli principle. Localization, in turn, 
reduces the overlap between the electron wave functions 
and, hence, the exchange interaction (a strongly localized 
system is believed to have a nearest neighbor antifer- 
romagnetic order [1^). This scenario also explains the 
large positive ma gnet oresistance observed in the vicin- 
ity of the ric HE 113 ■ Magnetic field aligns the spins, 
and again by the Pauli principle, drives the system to- 
wards the insulating phase. The localization transition 
at higher fields is, hence, shifted to higher densities. The 
proposed scenario also highlights the similarity between 
the localization transition in high mobility 2DEG and the 
Mott transition [j^. 

Note, that in contrast to all expectations, x the 
metallic phase depends on temperature down to 0.1 K. 
This dependence indicates the existence of a relevant 
energy scale considerably smaller than the Zeeman one 
(dashed vertical line in Fig. 0)). Such an energy scale 
may originate from localized spins which interact very 
weakly with each other. Quantification of the number of 
localized spins and their contribution to M [29j requires 
further study. 

Fig. El depicts the magnetization vs. magnetic field 
for various densities. The magnetization at densities 
close to Uc increases strongly with decreasing temper- 
ature. The magnetization is nonlinear, implying that 



To the best of our knowledge, with the exception of 
Refs. H [13 there is no reported experimental obser- 
vation of Stoner instability in 2D systems. In particu- 
lar, recent susceptibility measurements based on Sh-dH 
data 0, 0| , carried out down to the critical density of 
a superb sample {ric = 8 x 10^" cm~^), find a finite sus- 
ceptibility in the whole density range, in agreement with 
our result. We therefore turn to careful examination of 
the arguments used in 0, to claim the observation 
of such an instability. Both references rely on the mag- 
netoresistance measured as a function of in-plane mag- 
netic field. At high densities the resistance grows approx- 
imately quadratically with the field up to some density 
dependent field. Then it saturates or at least becomes 
weakly field dependent. It is believed, that at these den- 
sities the saturation field corresponds to full spin polar- 
ization. This large positive magnetoresistance is generic 
to all samples that show the so-called metallic phase in 
2D and is, hence, very likely to provide an important clue 
for the understanding of the latter phenomenon. The au- 
thors ofESEl have noticed that normalized magnetore- 
sistance curves, p{n, B)/ p{n,0) (magnetoconductance in 
the case of 0), measured at different densities, can be 
collapsed onto a single curve if the field is scaled by a 
density dependent field Bc{n) (we use the notation of lid 
Ref. Ilil utilizes somewhat different analysis in the same 
spirit). Moreover, for high densities, where magnetore- 
sistance saturation is observed, Bc{n) can be set to the 
saturation field. At lower densities (still above Uc) the 
saturation can no longer be observed but a scaling field, 
Bc{n), can still be found, so that the curves collapse one 
on top of another. The authors of these references no- 
ticed that Be vanishes approximately linearly when the 
density approaches ric, namely, Bc{n) oc n — ric- They 
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FIG. 6: Normalized magnetoresistance at different densities 
plotted versus scaled magnetic field, B/Bc{n). The field 
Bc{n), used to scale the data, is shown by dots in the inset. 
The extrapolation of the scaling field to zero at a finite den- 
sity, no, was used in to claim a ferromagnetic instability 
at n = no. 



Yet, at densities in the vicinity of the critical one we ob- 
serve almost free-spin like susceptibility, indicating nearly 
perfect compensation of the kinetic energy by the ferro- 
magnetic interaction. The possible relation between the 
large spin susceptibility at the critical density and the 
transition to strong localization calls for further theoret- 
ical and experimental studies. Understanding the role 
and nature of the localized spins might also turn to be 
important. 
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then argued that, since Bc{n) corresponds to full spin 
polarization at high densities, it should also correspond 
to full polarization at lower densities, where magnetore- 
sistance saturation is no longer observable. Following 
that logic all the way to the critical density, they con- 
cluded that the vanishing of Bc{n) at some finite density 
must indicate spontaneous polarization at zero field, i.e. 
the long awaited Stoner instability. We can not exclude 
Stoner instability in the superb samples used in 113, IitI 
but we can prove that the procedure used to conclude the 
instability is wrong. To that end we show in Fig. that 
our data obey the same scaling as inllQ. In anticipation 
of the same dependence of Bc(n) upon density as inllfilwe 
surmise Bc{n) oc n — riQ (inset to Fig. 6) and find that for 
uq — 1.15 X 10^^ cm^^ all our scaled magnetoresistance 
curves collapse onto a single curve (Fig. EJ. Following 
the arguments in we could have concluded Stoner 

instability at no, but our direct magnetization measure- 
ments at that density show finite susceptibility. The same 
fact is also refiected in Fig.UJs. If the high field magneti- 
zation is extrapolated to zero field (dashed line) one may 
have erroneously predicted instability. Carrying the mea- 
surements to smaller fields exclude that possibility. The 
wrong assumption of Refs. 0, is the identification of 
Bc{n) with a full polarization at all densities. Some of 
the authors of 17 later restricted their conclusion to the 
non-existing case of perfectly clean samples [3(1 |. 

In summary. Using a novel technique we were able 
to measure the thermodynamic spin magnetization of 
strongly correlated 2D electrons in a single 2D layer. Al- 
beit the substantial enhancement of the low temperature 
susceptibility, no ferromagnetic instability was observed. 
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